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ABSTRACT 

Risk and cost trade-offs have been simulated using a probabilistic method. The probabilistic method 
accounts for all naturally-occurring uncertainties including those in constituent material properties, 
fabrication variables, structure geometry and loading conditions. The probability density function of 
first buckling load for a set of uncertain variables is computed. The probabilistic sensitivity factors of 
uncertain variables to the first buckling load is calculated. The reliability-based cost for a composite 
fuselage panel is defined and minimized with respect to requisite design parameters. The optimization 
is achieved by solving a system of nonlinear algebraic equations whose coefficients are functions of 
probabilistic sensitivity factors. With optimum design parameters such as the mean and coefficient of 
variation (representing range of scatter) of uncertain variables, the most efficient and economical 
manufacturing procedure can be selected. In this paper, optimum values of the requisite design 
parameters for a predetermined cost due to failure occurrence are computationally determined. The 
results for the fuselage panel analysis show that the higher the cost due to failure occurrence, the 
smaller the optimum coefficient of variation of fiber modulus (design parameter) in longitudinal 
direction. 


2. INTRODUCTION 

Aerospace structures are complex assemblages of structural components that operate under severe and 
often uncertain service environments. They require durability, high reliability, light weight, high 
performance, at an affordable cost. Composite materials are potential candidates for meeting these 
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requirements. Composite materials possess outstanding mechanical properties with excellent fatigue 
strength and corrosion resistance. Their mechanical properties are derived from a wide variety of 
variables such as constituent material properties and laminate characteristics (fiber and void volume 
ratios, ply orientation, and ply thickness). These parameters are known to be uncertain in nature. 

In order to account for various uncertainties and to satisfy design requirements, knockdown (safety) 
factors are used extensively. These knockdown factors significantly reduce the design load of 
composite structures which result in substantial weight increase, but without a quantifiable measure of 
their reliability. This paper describes an alternate method which determines the structural reliability. 
This method is embedded in the computer code IP ACS (Integrated Probabilistic Assessment of 
Composite Structures) [ref. 1] for a comprehensive probabilistic assessment of composite structures. 
The schematic of IP ACS is shown in Figure (1). 

Since cost is a major driver for a structural design, optimization techniques should be sought to achieve 
the balance between maximum reliability and minimum cost. In this paper, reliability-based cost 
optimization is conducted to assess the risk and cost trade-ofls. 

3. IPACS COMPUTER CODE FOR RELIABILITY CALCULATION 

The IPACS Computer Code [ref. 1] has evolved from extensive research activities at NASA Lewis 
Research Center to integrate probabilistic structural analysis methods [ref. 2] and computational 
composite mechanics [ref 3], The composite micromechanics, macromechanics and laminate theory 
are embodied in ICAN [ref 3], IPACS consists of two stand-alone computer modules: PICAN and 
NESSUS. PICAN simulates probabilistic composite mechanics [ref. 4], NESSUS uses the 
information from PICAN to simulate probabilistic structural responses [ref. 5], Direct coupling of 
these two modules makes it possible to simulate uncertainties in all inherent scales of the composite - 
from constituent materials to the composite structure including its boundary and loading conditions as 
well as environmental effects. It is worth noting that special reliability algorithm FPI (Fast Probability 
Integrator) (ref. 6) is used instead of the conventional Monte Carlo simulation [ref. 7], to achieve 
substantial computational efficiencies which are acceptable for practical applications. Therefore, 
probabilistic composite structural analysis becomes feasible which can not be done traditionally, 
especially for composite structures which have a large number of uncertain variables. The results from 
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IPACS analysis indude probability distribution function of structural response, reliability for a design 
criterion, and the probabilistic sensitivity factors of the primitive variables to the structural response and 
structural reliability. 

4. PROBABILISTIC SENSITIVITY INFORMATION 
FOR OPTIMIZATION 

The commonly used sensitivity in a deterministic analysis is the performance sensitivity, dZJdXi, which 
measures the change in the performance Z due to the change in a design parameter X,. This concept is 
extended to probabilistic analysis to define the probabilistic sensitivity which measures the change in 
reliability relative to the change in each random variable. The failure probability for a given 
performance is defined in equation 1 (ref. 8). 

P f = W 1 


where fi is the r eliab ility index; O is the cumulative distribution function of a normally distributed 
random variable. Probabilistic sensitivity factor (SF;) for i* random variable is defined in equation 2. 
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where u is the most probable failure point of a limit state function in a unit normal probability space. 
The sensitivity of design parameters to structural reliability is another useful information for controlling 
and adjusting design parameters from manufacturing to obtain the "best" benefit with minimum 
al tpr^tirm The sensitivity of the reliability to the mean of a normally distributed random variable X; 
can be represented by equation 3 (ref. 8). 



where m t and cr, are the mean and standard deviation of random variable X; respectively; Similarly, 
the sensitivity of the reliability to the standard deviation of a normally distributed random variable Xj 
can be computed from equation 4. 

dfi _ SFj u’j = u? ^ 

0 (j j cy j ex * 

With this information, optimization with reliability considerations can be achieved as will be 
demonstrated later. 

5. MINIMIZATION OF RELIABILITY-BASED COST 

A major design goal is to achieve the balance between maximum reliability and minimum cost. The 
criterion which addresses both reliable component and cost simultaneously is the cost function shown 
in equation S. 

Ct = Cj P f Cf 5 


where Ct is total cost; Ci is manufacturing cost; P f is failure probability; Cf is the cost incurred 
due to structural Mure. Manufacturing cost can be represented by equation 6. 

Ci = Z Cj(Pj) + Co 6 

i-i 


where N is the number of distribution parameters; C> ( Pj) is manufacturing cost for j* distribution 

parameter, Co is other unrelated constant cost (assume zero in this paper). Total cost can be 
minimized when 


dCr_ 

~dPj 


= 0 


j- 1, N 
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Therefore, optimization can be achieved by solving a system of nonlinear algebraic equations as shown 
in equation 8. 


dCjfPj) dty-Pl dp_ = 

dPj dp d Pj 


For a normally distributed random variable, dp / dp j can be calculated by equations 3 and 4. 


6. DEMONSTRATION FOR OPTIMIZATION OF COMPOSITE 
FUSELAGE PANEL STRUCTURES 

Reliability-based cost minimization is demonstrated using a lower side panel of a composite fuselage 
structure shown in Figure (2). The fuselage panel is 60 in. x 90 in. with a radius of 122 in. There are 
four longitudinal stringers, 1.49 in. x 0.304 in., four J-frames, 5.10 in. x 0.140 in., and a 0.096 in. thick 
drin with 0. 191 in. thick pads under the stringers and 0. 124 in. thick pads under the frames. The width 
of the pad for stringer and frame are 3.57 in. and 2.85 in. respectively. The panel is made of 
graphite/epoxy composite. The constituent material properties, their assumed probabilistic distribution 
and coefficient of variations (representing range of the scatter) are summarized in Table 1. The 
corresponding fabrication variables used to make the composite panel are summarized in Table 2. 
Although only normally distributed random variables are considered in this example, the computer 
code IP ACS can handle random variables with other distribution types such Lognormal or Weibull. 
The reliability of the structure is evaluated for design criterion that requires the buckling load be greater 
than 30% of the ultimate design load. 

The probability density function of the buckling load is shown in Figure 3. Based on a reliability of 
0.999, the ultimate design load will be 5000 lb/in. (1500 Ib/in. divided by 0.30). The probabilistic 
sensitivity factors of the uncertain variables to 0.001 probability of the buckling load are listed in Figure 
4. It is found that fiber volume ratio has the highest influence on the scatter of the buckling load, 
followed by the skin thickness, pad thickness under stringer and fiber modulus in longitudinal direction. 
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The probability of Mure for different coefficients of variation of Em with the rest of distribution 
parameters unchanged is shown in Figure 5. It is noticed that probability of Mure increa<ses while 
coefficient of variation increases. It means that lesser quality product has higher probability of Mure. 

To minimize the total cost in equation 5, initial cost Cj is defined to be 

Ci = A( Pl ) * p, + B(p 3 ) * D * (1-pJ 9 


where p l and p 2 are the coefficients of variation of fiber modulus (Em) and matrix modulus (Em) 
respectively; p 3 represents the mean value of fiber volume ratio. A( p,) is the cost for fiber defined 
by equation 10. 

A(pJ = Ao * (1.25 -p]) 10 


where Ao is set to be $35/lb to represent typical cost for graphite fiber. B( p,) in equation 9 is the 
cost for epoxy matrix defined by equation 1 1 . 

B(P,) = Bo * (1-25 -p 2 2 ) 11 


where B 0 is set to be $25/Ib to represent typical cost for epoxy matrix. D in equation 9 is the ratio of 
the weight densities between epoxy matrix and graphite fiber. The cost function defined by equation 9 
only accounts for the cost to manufacture fiber and matrix material. Other manufacturing cost is 
assumed to be constant and is set to be zero in this assessment. To simplify the optimization procedure 
and for demonstration purposes, only the coefficient of variation of fiber modulus in longitudinal 
direction is chosen to be the design parameter. Substituting equation 9 into equation 8, optimization 
can be achieved by solving the following equation. 
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&b(-p) dp 
dp dp, 


Letting Cf be $ 15000/lb, the manufacturing cost Q and the weighted failure cost Pf*C F at different 
coefficient of variation (COV) for Em are plotted in Figure 6. It is noticed that Q decreases while 
COV of Em increases because low quality fiber (wider scatter of Em) costs less to manufacture. On 
the other hand, the weighted failure cost increases since more weight for the failure cost is considered. 
Optimum pi (COV of Em) is calculated using equation 12. The minimized reliability-based total cost is 
calculated by equation 5. In Figure 7, total cost is computed at different COV of Em. It is seen that at 
optimum COV of Em (0.05), the manufacturing cost and weighted cost due to Mure is balanced. 
However, the increase in the weighted cost due to Mure is greater than the reduction in the 
manufacturing cost when COV of Em is greater than 0.05. The optimum COV of Em for different 
cost due to failure occurrence is computed as shown in Figure 8. It is seen that the larger the cost due 
to Mure, the smaller the optimum COV of Em. 

From this assessment, it demonstrates that selection among posable arrangements for optimum cost 
with reliability consideration can be achieved using the probabilistic method adopted in the computer 
code IP ACS. 


7. CONCLUSIONS 

Risk and cost trade-offs for polymer composite structures are assessed and demonstrated. The 
assessment is accomplished by integrating optimization of reliability-based cost concept into the IP ACS 
(Integrated Probabilistic Assessment of Composite Structures) computer code. The optimization is 
achieved by solving a system of nonlinear algebraic equations whose coefficients are functions of 
probabilistic sensitivity factors. Optimum design parameters such as the mean and coefficient of 
variation (representing range of scatter) of uncertain variables are essential information for the selection 
of the most efficient and economical manufacturing procedure. In this paper, the reliability-based cost 
for a composite fuselage panel is minimized with coefficient of variation of fiber modulus in 
lon gitudinal direction as the design parameter. The minimum total cost for a 15000 Mure cost is 
found when COV of Em is equal to 0.05. Therefore, the balance between the maximum reliability and 
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mini mum cost is achieved. The optimum COV of Em for different failure cost is also computed. It is 
found that the larger the failure cost, the smaller the optimum COV of Em- 
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9. SYMBOLS 


Em: fiber modulus in longitudinal direction 

Ec 2: fiber modulus in transverse direction 

G fl2 : in-plane fiber shear modulus 

Gf23: out-of-plane fiber shear modulus 

Vn2: in-plane fiber Poisson's ratio 

V123: out-of-plane fiber Poisson's ratio 

E,,,: matrix elastic modulus 

Gn,: matrix shear modulus 

v m : matrix Poisson's ratio 

fvr: fiber volume ratio 

pdf: probability density fimction 

<I> : cdf of Normal Random Variable 

Z : performance function 

fi : Reliability Index 

P f : Probability of Failure 

SF : Safety Factor 

u : Most Probable Point 

N : Number of Distribution Parameters 

Co ■ Constant Cost 

Pj : j* Distribution Parameter 

Cr : Reliability-based Total Cost 

Ci : Manufacturing Cost 

C F ■ Cost due to Failure 

Xi : i th random variable 

m, : mean of i th random variable Failure 

cov: coefficient of variation 

pi: COV ofEfii 

P2: COVofEn, 

p3 _ mean value of fiber volume ratio 

Ao'. typical manufacturing cost for graphite fiber 

Bo: typical manufacturing cost for epoxy matrix 

D : ratio of weight densities between matrix and fiber 

a i : standard deviation ofith random variable 

Cj (Pj) '■ Manufacturing Cost for j 4 Distribution Parameter 
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Table 1, The Statistics of Fiber and Matrix Properties for Graphite-Epoxy Composite 


PROPERTY 

SYMBOLS 

UNIT 

DISTRIBUTION TYPE 

MEAN 

COEFFICIENT OF 
VARIATION 
(percentage of scatter) 

Eni 

Msi 

Normal 

31.0 

0.05 

Ee, 

Msi 

constant 

2.0 

0.00 

Gfl2 

Msi 

constant 

2.0 

0.00 

Gq3 

Msi 

constant 

1.0 

0.00 

V i2 

— 

constant 

0.2 

0.00 

Vj3 

— 

constant 

0.25 

0.00 

Em 

Msi 

Normal 

0.5 

0.05 

G m 

Msi 

constant 

0.185 

0.00 

Vm 

— 

constant 

0.35 

0.00 


Table 2. The Statistics of Fabrication Variables 


PROPERTY SYMBOLS 

UNIT 

DISTRIBUTION 

TYPE 

MEAN 

COEFFICIENT OF 
VARIATION 
(percentage of scatter) 

Fiber Volume Ratio 

— 

Normal 

0.600 

0.05 

Void Volume Ratio 

- 

constant 

0.020 

0.00 

Ply Misalignment Angle 

deg. 

constant 

0.000 

0.00 

Skin Thickness 

in 

Normal 

0.096 

0.04 

Stringer Thickness 

in 

Normal 

0.304 

0.04 

Frame Thickness 

in 

Normal 

0.140 

0.04 

Pad Thickness (Stringer) 

in 

Normal 

0.191 

0.04 

Pad Thickness (Frame) 

in 

Normal 

0.124 

0.04 
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Figure 2. - Lower Side Panel Component 
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Figure 5. - Probability of Failure as a Function of COV of Efl 1 



Figure 6. - Normalized Partial Cost as a Function of COV of Efl 1 
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Cj : Manufacturing Cost 
Pf : Failure Probability 
Cp Failure Cost 
C T : Total Cost 
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Figure 7. - Normalized Total Cost as a function of COV of Efl 1 
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Figure 8. - Optimum COV of Fiber Modulus Efl 1 as a Function 
of Normalized Cost due to Failure Occurrence 
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